Getting the most from the statistical analysis of solar neutrino oscillations 
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We present a thorough analysis of the current solar neutrino data, in the context of two-flavor 
active neutrino oscillations. We aim at performing an accurate and exhaustive statistical treat- 
ment of both the input and the output information. Concerning the input information, we analyze 
81 observables, including the total event rate from the chlorine experiment, the total gallium event 
rate and its winter-summer difference, the 44 bins of the Super-Kamiokande (SK) energy-nadir elec- 
tron spectrum, and the 34 day-night energy spectrum bins from the Sudbury Neutrino Observatory 
(SNO) experiment. We carefully evaluate and propagate the effects of 31 correlated systematic 
uncertainties, including 12 standard solar model (SSM) input errors, the 8 B neutrino energy spec- 
trum uncertainty, as well as 11 and 7 systematics in SK and SNO, respectively. Concerning the 
output information, we express the x 2 analysis results in terms of "pulls," embedding the single 
contributions to the total x 2 coming from both the observables and the systematics. It is shown 
that the pull method, as compared to the (numerically equivalent) covariance matrix approach, is 
not only simpler and more advantageous, but also includes useful indications about the preferred 
variations of the neutrino fluxes with respect to their SSM predictions. Our final results confirm 
the current best-fit solution at large mixing angle (LMA), but also allow, with acceptable statisti- 
cal significance, other solutions in the low-mass (LOW) or in the quasi-vacuum oscillation (QVO) 
regime. Concerning the LMA solution, our analysis provides conservative bounds on the oscilla- 
tion parameters, and shows that the contribution of correlated systematics to the total \ 2 is rather 
modest. In addition, within the LMA solution, the allowed variations from SSM neutrino fluxes are 
presented in detail. Concerning the LOW and QVO solutions, the analysis of the pull distributions 
clearly shows that they are still statistically acceptable, while the small mixing angle (SMA) solution 
could be recovered only by ad hoc "recalibrations" of several SSM and experimental systematics. A 
series of appendices elucidate various topics related to the x 2 statistics, the winter-summer differ- 
ence in GALLEX/GNO, the treatment of the SK and SNO spectra, and a quasi-model-independent 
comparison of the SK and SNO total rates. 

PACS numbers: 26.65.+t, 13.15.+g, 14.60.Pq, 91.35.-x 



I. INTRODUCTION 

The data from the Homestake §, SAGE Jp, GALLEX/GNO |,|], Kamiokande [|, Super-Kamiokande (SK) [§, 
and Sudbury Neutrino Observatory (SNO) (t]7||, Pj experiments have consistently established that electron neutrinos 
emitted from the Sun JTo| undergo flavor transitions to the other active states (y^ or v T ). Neutrino oscillations [pd| , 
possibly affected by matter effects in the Sun or in the Earth represent a beautiful explanation of such transitions. 

Assuming the simplest scenario of two-family oscillations among active neutrinos, an important task for the next 
future is to refine the current constraints on the neutrino squared mass difference 5m 2 = m| — mf > and on the 
mixing angle Orz € [0,7r/2]. In order to accomplish this task, one needs: (a) new or more precise measurements; (b) 
accurate calculations of the v e survival probability P ee (Sm 2 , 612) and of related observable quantities; and (c) powerful 
statistical analyses to compare the (increasingly large) solar v data set with theoretical expectations. 

The point (a), not discussed in this work, will soon be addressed by the decisive reactor v e experiment KamLAND 
[13 1, as well as by the solar v experiments which are currently running ^, 0], being restored ||, or in construction 
[14|. Concerning the point (b), since the current numerical and analytical understanding of the oscillation probability 
(and related observables) is quite mature in the whole (5m 2 , 612) plane, we will only make a few remarks when needed. 
In this paper, we rather focus on point (c), aiming to an exhaustive statistical analysis including all known observables 
and uncertainties in input, and providing very detailed information in output, in order to better appreciate the current 
status of the solutions to the solar v problem in terms of active flavor oscillations. Although some details will be 
specific of solar v data, the analysis method that we discuss is quite general, and can be easily extended to any kind 
of global fit. 

The structure of our paper is the following. In Section II we discuss the equivalence between the \ 2 approaches in 
terms of the covariance matrix and of "pulls" of observables and systematics. In Sec. Ill we describe the input and 
output of the pull approach, as applied to the analysis of 81 solar neutrino observables and of 31 input systematics, in 
the context of 2v active oscillations. In Sec. IV we discuss the x 2 analysis results in terms of multiple allowed regions 
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in the mass-mixing parameter space (Sm 2 , tan 2 6*12) and in terms of the associated pull distributions, which provide 
additional information about the relative likelihood of the various solutions and about the allowed deviations from 
standard solar model (SSM) fluxes. We draw our conclusions in Sec. V. More technical (but sometimes substantial) 
issues are discussed in a series of Appendices, which deal with the x 2 statistics (App. |A|), the winter-summer asymmetry 
in GALLEX/GNO (App. |b|), the treatment of the SK spectrum uncertainties (App. |Cj), the analysis of the SNO data 
(App. |j5|), and a quasi- model-independent comparison of SK and SNO total rates (App. ||). 

As a conclusion to this Introduction, we would like to stress that deepening the statistical analysis and improving 
the evaluation of the uncertainties is an important task in neutrino oscillation physics, just as it happens (or happened) 
in other areas of "precision" physics. Indeed, after the observation of two large oscillation effects (the disappearance 
of atmospheric za, and of solar v e , and their upcoming tests at long-baseline accelerator and reactor experiments), we 
are likely to face an era of delicate searches for smaller effects related, e.g., to the angle #13, to leptonic CP violation, 
or to subleading contributions induced by nonstandard v states or interactions. Moreover, one should not forget that, 
so far, there is no direct evidence for a vacuum oscillation pattern (disappearance and reappearance of a specific 
flavour) or for matter effects in the Sun or the Earth. Such effects might well generate only small signals in present or 
planned experiments, and any effort should be made in order to quantify them (if any) with accurate analyses. From 
this viewpoint, we think that our thorough analysis can add valuable information and useful technical tools to other 
solar v fits 0, |, ||, ||, @ HI 13 that appeared soon after the release of the SNO neutral current data flSfl . 



II. TWO EQUIVALENT WAYS OF DEFINING THE x 2 FUNCTION 

Let us consider a set of N observables {Rn}n=i,...,N with their associated sets of experimental observations {R°* pt } 
and theoretical predictions {R^ c °}. In general, one wants to build a \ 2 function which measures the differences 
i?° xpt — i?J; hcor in units of the total (experimental and theoretical) uncertainties. This task is completely determined 
if, for any difference i?° xpt — i?^ hoor , one can estimate an uncorrelated error u n , and a set of K correlated systematic 
errors c„ induced by K independent sources, namely 

R^ t -R t ^±u n ±c 1 n ±cl...±c« (n=l,...,N) , (1) 

with 

p(u n ,u m ) = S nm , (2) 
p(ci> c t,) = s kh V(n,m) , (3) 

where p represents the correlation index. 1 

Given the input numbers in Eq. (|l|), two \ 2 definitions have been basically used i n gl obal analyses (with some 



varia nts o r combinations). They will be referred to as the "covariance" approach (Sec. II A) and the "pull" approach 



(Sec. [IB). Although seemingly different, the two approaches are strictly equivalent (Sec. II C). The pull approach 



however, proves to be much more advantageous, and will be used throughout this paper. 

A. The covariance approach 

In the "covariance approach," one builds the (covariance) matrix of squared errors as 

K 



a nm ~ $nm,Un u m + ^ ] C n C m > (^) 



k=l 



then inverts it, and evaluates the quadratic form 

N 



2 _ \ * /poxpt r>thcor\r 2 1 — 1 f jjexpt ptheor\ (e\ 

Acovar / , l^Vi ■ n "n )l u nm\ (■ n "m ) • \° I 



1 The error represents the shift of the n-th observable induced by a +1<t variation in the fc-th systematic error source. Linear 
propagation of errors is assumed, namely, possible ±1ct asymmetries and second-order systematic effects oc c*cjj are consistently 
neglected in computing the uncertainties of the difference _R^ xpt — i?* cor . 
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This approach, proposed in |2l| for the data available at that time, has been later used in the majority of solar v 
analyses of total events rates, with some variants related to the treatment of the 8 B v flux (free or SSM) and to the 
separation of spectral and total rate information in the SK data. 



B. The pull approach 

The alternative "pull approach" embeds the effect of each independent fc-th source of systematics through a shift 
of the difference i?° xpt — i?^ heor by an amount — £fcc£, where £/. is a univariate gaussian random variable 2 

K 

(iC Pt - ^ h0 ° r ) "> (iC Pt - ^ hC ° r ) - E ■ ( fi ) 

fe=l 

The normalization condition for the £fc's is implemented through quadratic penalties in the x 2 , which is then minimized 
with respect to all £fc's, 



Xpull = min 



N ( ncxpt nthcor C, Jc\ K 

n=l \ / fc=l 



(7) 



Denoting as £ fc ("pulls" of the systematics) the values of the £fc's at the minimum, and defining the "pulls" x n of the 
observables as 

_ _ J% xpt -« hror + Ef=i? fc <£) 

the value of Xp U n ' s tben split into two diagonalized pieces, embedding the contribution to the \ 2 from the residuals 
of the observables and of the systematics, 

.2 _ ,,2 , ,,2 



Xpull Xobs Xsys (9) 
N K 

= D£+E£- (10) 



1 —2 

fc=l 



The pull approach has often been used by the SK collaboration in their day-night spectral analysis, in combination 
with the covariance method for non-SK data [p2|. More recently, the SK energy-nadir spectrum has been analyzed 
through a mixed pull-covariance approach [E3T The link between the covariance and pull method, discussed in the 
next section, is also mentioned in passing in Rcf. p4| . To our knowledge, however, a complete analysis of solar neutrino 
data in terms of Xp U n nas n °t been performed, prior to the present work. 



C. Comparison and equivalence of the covariance and pull approaches 

It is perhaps not generally known that, although seemingly different, the covariance and pull approaches are strictly 
equivalent, 

X 2 = X 2 u (11) 

Acovar — Apull 5 \ J 

Our proof of the above identity is given in Appendix [A[ 3 Given the equivalence in Eq. (|ll]) , the choice between the 
covariance and the pull approach must be dictated by their relative merits. 



is conventional. It amounts to attribute all the shifts to the theoretical estimate: i?,^ llcor — * 

nthcor i c k 

^ fc=1 \ — k o o o 

We have recently realized that Eq. ( |ll[) an d its implications have been also discussed [123] and are routinely used K6l E7|] in the context 
of parton density distribution fitting |28|. Closely related results have also been recently found in the context of cosmic microwave 
background data fitting [E9[ . Although the connection between covariance matrix and pulls appears thus to be an ubiquitous result in 
physics data analysis, we have been unable to trace explicit references to Eq. (|ll]) prior to Ref. 
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In particle physics, the covariance approach is typically used either when the experimental collaborations provide 
detailed information about the correlation matrix (as in the case of the LEP Electroweak Working Group |3(| ) , or when 
N <C K (as in the case of solar neutrino fits to total rates only ^J). However, for increasing N the approach becomes 
increasingly complicated. The inversion of large N x N covariance matrices, besides being numerically tricky, can 
make it difficult to fully understand the results of global analyses. Indeed, the current solar or atmospheric neutrino 
data fits, involving N ~ O(10 2 ), are getting close to their manageability limits in terms of covariance matrices. The 
situation might become even more problematic in future high-statistics experiments, such as the neutrino factories 
from muon storage rings or superbeams, where the oscillation parameters will be inferred from the analysis and 
comparison of densely binned and correlated (anti) neutrino event spectra. 4 

The pull approach is clearly more practical than the covariance one when K -C N . In fact, the minimization in 
Eq. (^) leads to a set of K linear equations in the £ fe 's, and to an associated K x K matrix inversion, rather than the 
N x N covariance matrix inversion (see App. [A]). Moreover, the final decomposition in terms of pulls of observables and 
systematics [Eqs. ([)]) and (1Q|)] allows to trace the individual contributions to the x 2 , and to easily detect anomalously 
large residuals. 5 In general, this method is useful to gauge the mutual agreement of data in a global fit, or to diagnose 
tension among data (if any), for any given point in the model parameter space [(6m 2 , tan 2 #12) in our case]. The 
analysis of the Xohs an< ^ Xs ys components [Eq. (||)] can also be useful to trace possible sources of good or bad fits. 

Given the advantages of Xp U \i m cases where K <C N, we have redesigned the statistical analysis of solar neutrino 
oscillations in terms of pulls, and applied it to the current data set (where N = 81 and K = 31), as discussed in the 
next section. 



III. INPUT AND OUTPUT FOR THE Xp U n ANALYSIS 

In this section we describe the main input and output quantities, related to the Xp U n analysis of solar neutrino data. 
In input we consider a set of N = 81 observables R n (with associated uncorrelated errors u n ), and a set of K = 31 
sources of correlated systematic errors in part related to the SSM and in part to the experiments. 6 In output we 
consider the total Xpuii; its decomposition in individual pulls, and the shifts of the neutrino fluxes from their SSM 
value. 



A. Input observables and uncorrelated errors 

The first two observables in our list are the Chlorine total rate [Q, 

= 2.56 ± 0.23 SNU , (12) 
and the average Gallium total rate (SAGE g + GALLEX/GNO @), 

jjjexpt = 7Q g ± 4 4 SNU ^ ( 13 ) 

where the errors include the statistical and experimental systematic contributions to the uncorrelated errors uc\ 
and «Qa, respectively. In the analysis, the corresponding cross section uncertainties must also added in quadrature. 
Following the suggestion in ||]|, the cross section error components ARx.i for X =(C1, Ga) are first added linearly and 
then quadratically into low (L) and "high" (H) energy parts [L =(pp, pep, Be, N, O) and H =(B, hep), respectively] 7 



\i£L 




u\ (cross section) = £ Ai#~ ] + ( A R%? ) . (14) 



4 The pull approach has been recently applied to prospective studies in this context |3l| . 

5 The usefulness of the pull distribution diagnostics has been already recognized in many areas of physics, including solar neutrino physics 
Jit] , [32| and electroweak precision physics ]jo| . However, in Refs. Jis| , |30| p^ ] the pulls are denned without the shifts in Eq. (g), and are 
thus all correlated to each other. 

6 Systematic error sources are defined as "correlated" if they act upon two or more observables at the same time. Systematics which act 
upon one observable only (e.g., the v-Cl absorption cross section error) simply contribute quadratically to the uncorrelated error for 
that observable (the CI t ota l rate, in the example). 

7 The prescription in Eq. (J14) is intermediate between the extreme possibilities of quadratic sum Ell and linear sum ]34[ over all flux 
components. We observe that this prescription is not only justified by the physics of the u e absorption processes in CI and Ga [B3| , but 
also by the effective separation of the CI and the Ga solar v response functions into two "L" and "if clusters in the energy aomain 
§§§. 
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For X=G\, the cross section error components are evaluated as ARx.i = Rx,i A In Cx.i, where Rx.i are the (oscillated) 
rate components, and the fractional la cross section uncertainties AhxCxi can be taken from the compilation in 
p7j . For A=Ga, the value of ARx,i is computed by taking 1/3 of the variations induced by the ±3cr perturbed cross 
sections |j8| on the i-th Ga rate component for each point of the oscillation parameter space, as suggested in f^Sf. 8 

Our third observable is the winter-summer (W — S) rate difference |3!| measured in GALLEX/GNO 0], here 
introduced for the first time in the oscillation analysis. This datum is described in detail in App. |b|. 

The SK experiments provides 44 observables, in terms of (binned) absolute event rates for the energy-nadir differ- 
ential spectrum of electrons Our treatment of the SK spectral information is described in detail in Sec. |c[ 

The set of solar v observables is completed by the 34 day-night energy spectrum bins from the SNO experiment 
||, ||, which include contributions from v elastic scattering (ES), charged current (CC) and neutral current (NC) 
interactions, and backgrounds [Eol. Our treatment of the SNO spectral information is described in App. O. 



B. Input correlated systematics 

The 31 sources of correlated systematics include 12 uncertainties related to SSM input, the 8 B v energy shape 
uncertainty, 11 SK error sources and 7 SNO error sources. 

Concerning the SSM input, we take from the central values for the fluxes (^pp, 'I'pep, $Be, 3>b, $n, $o), but 
rescale $hep from the value in 0] (9.3 x 10 3 cm" 2 s" 1 ) to our default value 

$ hcp = 8.3 x 10 3 cm" 2 s _1 , (15) 

according to the recent evaluation of the associated Shep factor in jl2| . 

The SSM also embeds a set of eleven sources X k of correlated uncertainties (the cross section factors Su, S13, 
S34, Si, 14, 517, the Be capture cross section Cb c , the Sun luminosity, metallicity Z/X, age, opacity, and element 
diffusion), with fractional uncertainties AhiXk, as listed in |37j |. With respect to the compilation in j37j], we update 
A In Z/X = 0.061 from jyj , and we add X12 — S'hcp, with A In Shr., n = 0.3. 9 The effects of such sources of uncertainties 
on the neutrino fluxes <&i are characterized by log-derivatives [|10| ^H , a,t = <91n$i/<91n Afc, as compiled in [j37| . 
Concerning X12, the only nonzero log-derivative is «hep,i2 = 1- 

The collective effect of the SSM sources of systematics X k amounts to shift the neutrino fluxes as 

^^<P t (^ + J2^ k a lk A\nX^j , (16) 

where the £ fc 's, penalized by the quadratic term £ic m the expression of Xpuii' are minimized away in the fit. 

Notice that the above equation is the linearized form of the power laws connecting each flux to the AVs [[Io| ^3| . 
Such linear form satisfies the luminosity constraint for the fluxes fl44| by construction, due to the sum rule discussed 
in Ref. f| |t| (see also g|). 

Within the pull approach, the shifts in Eq. (|l6|) are easily propagated to all theoretical predictions i?^ heo . In 
particular, if i?*j h °° is the z-th flux contribution to i?^ hco , then the associated correlated shift from the fc-th source is 
£fc c n i = ikR t nT >a ikAhi X k . The net effect of the shifts in Eq. ( |T6| ) is thus the generation of correlated errors on the 
i?„'s, which is strictly equivalent to the construction of the astrophysical error matrix defined in the earlier covariance 
approach [^l], |37j , as also noticed in . 

The 13th source of correlated systematics in our list is the 8 B v spectrum shape uncertainty [fl6[| around the current 
"central" spectrum f47|] , which affects all the 81 observables 10 at the same time. In the absence of oscillations, we 
estimate that a +lcr perturbation of the 8 B v spectrum (in the direction of higher v energies) generates a +2.2% and a 
+1.7% increase of the $s component of the i?Ga and Rci theoretical rates, respectively (the evaluation is repeated in 
each point of the mass-mixing plane). The corresponding shifts for the SK and SNO spectra are evaluated in App. ^ 
and ID], respectively. 

There are also eleven sources of correlated systematics, which affect only the SK spectrum. They include the 
SK energy scale and resolution uncertainties, an overall SK rate offset, and eight sources of systematics separately 
affecting the eight energy bins, with full correlation in nadir || (see App. Q). 



We conservative v assume the largest between the +1ct and —la asymmetric Ga cross section errors. 
The authors of [12| quote an uncertainty of ~ 15% for S^ ep , that we conservatively double to 30%. 

With the possible exception of the Ga winter-summer difference, where its effects cancel to a large extent, and can be safely neglected 
as compared with the rather large statistical error (see App. ph. 
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Finally, there are seven sources of correlated systematics, which affect only the SNO spectrum. They include: the 
uncertainties affecting the SNO energy scale and resolution, the event vertex reconstruction, the neutron capture 
efficiency, the neutron and low-energy (LE) background estimates and the interaction cross sections. See App. |5] for 
more details. 

C. Output 

As output of the pull analysis, we get the function Xp ull (5m 2 , tan 2 6*12) (essential to identify absolute and local 
minima and to draw confidence level contours), as well as other useful statistical indicators. 

Concerning Xp U u> f° r an y fi xe d point in the parameter space {5m 2 , tan 2 #12), the goodness-of-fit test requires Xp U n ~ 
N [N = 81) for an acceptable fit. Further information can be gained by splitting Xpuii m t° t ne separate contributions 
X 2 bs and Xsys [Eqs. || and [To| ], obtained by summing up the squared pulls of the N = 81 observables \x n , see Eq. (]§)] 



and of the K = 31 systematics [£ fe , see Eq. (A6)]. The larger the value of Xsys> the more the fit tends to "stretch" 
one or more correlated systematics to get a better agreement between data and expectations. Apart from global 
features, the analysis of the pull sets {x n } and allows to quantify individual contributions to the x 2 > which, 

if anomalously large, might be indicative of problems either in the theoretical predictions or in the experimental 
measurements. Therefore, we think it useful to present, besides the global values of Xp U n = Xobs + Xsysi also some 
selected lists of pulls. 

Finally, it is useful to isolate the twelve SSM systematic pulls {£fc}fc=i,...,i2 which, on the basis of Eq. (|l^), allow 
to derive the induced neutrino flux shifts from the SSM central values, namely 

A$ 12 



= ^£ fc a lfc AlnX fe . (17) 



k=l 



These shifts provide valuable (and luminosity-constrained) information about the preferred departures from the SSM 
within the various oscillation solutions to the solar neutrino problem. 
Summarizing, we will show and discuss results about pulls, 



about x 2 values, 



{#n}n=i,...,8i — pulls of the observables , (18) 
{£,k}k=i 31 = pulls of the correlated systematics , (19) 



(20) 



Xobs 


81 




n=l 




31 


Xsys 


= ^2~tk > 




fe=l 




_ 2,2 


Xpull 


Xobs ' Xsys 



(21) 
(22) 

and about fractional shifts from the SSM predictions, 

A$i/$ j; = v flux shifts . (23) 

IV. RESULTS OF THE Xp U n ANALYSIS 

In this section we start by describing the global results of the Xp U n analysis, and then we break down such results 
at increasing levels of detail. 

A. Global results 



The global results of our solar v oscillation fit are summarized in Table | and in Fig. [j]. In Table | we report the 
(5m 2 , tan 2 6x2) coordinates of the best-fit point [so-called large mixing angle (LMA) solution] and of the three deepest 
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TABLE I: Positions and values of the absolute minimum (LMA) and of three relevant local minima (LOW, QVO, SMA) 
of Xpuii! together with the separate contributions from pulls of observables (Xobs) an d of correlated systematics (Xsys)- The 
corresponding A^ 2 variations are also given. 



Solution 


Sm 2 (eV 2 ) 


tan 2 #i2 


Xobs 


AxiL 


Xsys 


AXs 2 ys 


Xpull 


AXpuii 


LMA 


5.5 x 1(T 5 


0.42 


71.3 




2.1 




73.4 




LOW 


7.3 x 10" 8 


0.67 


79.7 


8.4 


4.1 


2.0 


83.8 


10.4 


QVO 


6.5 x KT 10 


1.33 


74.9 


3.6 


6.3 


4.2 


81.2 


7.8 


SMA 


5.2 x 1(T 6 


1.1 x 10~ 3 


83.1 


11.8 


13.8 


11.7 


96.9 


23.5 



(local) Xpuii m i n i ma m the regions of low Sm 2 (LOW), quasivacuum oscillations 11 (QVO), and small mixing angle 
(SMA). 

Concerning the goodness-of-fit test, we remind that, at the absolute minimum, one expects the total x 2 to be in 
the ±lcr range iV DF ± ^2N DF @. In our case (N DF = 81 - 2), it is Xp U n = 73.4 at the LMA best-fit point, well 
within the expected range 79 ± 12.6. Also the LOW and especially the QVO solutions have acceptable values of Xp U n> 
while the SMA value appears to be significantly larger than expected. The pull analysis will confirm that the LOW 
and QVO solutions are still viable, while the SMA solution is no longer statistically acceptable. Notice that at the 
LMA point, most of the contribution to Xp U n comes from pulls of observables (Xobs) rather than systematics (x 2 ys )- 
All the other solutions in Table | show an increase of both Xpuii and x 2 ys i implying an increasing departure of the 
theoretical predictions from the data and of the systematics offsets from zero. 

Concerning the parameter estimation test 12 (based on Ax 2 u u variations around the minimum) Figure 1 shows the 
results of our global analysis in the usual mass-mixing plane. The confidence level isolines are drawn at Ax 2 = 4.61, 
5.99, 9.21, and 11.83, corresponding to 90%, 95%, 99% and 99.73% joint probability regions for the two (Sm 2 , tan 2 6> 12 ) 
parameters. The QVO and LOW parameters are still acceptable at the 99% and 99.73% C.L., respectively, while the 
SMA parameters are basically ruled out. Our LMA bounds appear to be (sometimes significantly) more conservative 
than in other recent analyses §, |, || [l|, [TJ, [jj], || Q. 13 In particular, at the 99.73% C.L. we derive from 
Fig. [j] that: (i) maximal mixing is marginally allowed in the LMA region, and (ii) the highest Sm 2 allowed values 
hit the region independently disfavored by CHOOZ data |5!|. 14 We think that the detailed treatment of all known 
uncertainties (and of their propagation to all relevant experimental observables) plays a role in such different results, 
also for non-LMA solutions. Concerning the LOW solution, we note that the inclusion of the winter-summer datum 
from GALLEX/GNO contributes to decrease its likelihood in our analysis. 



B. Separating experimental bounds 

Figure|| show the decomposition of the global results into contributions from the CI experiment (total rate), from 
the Gallium experiments (total rate and winter-summer difference), from the SK energy- nadir spectrum (44 bins), 
and from the SNO day- night spectrum (34 bins). In each panel, the results are shown in terms of allowed regions, for 
the same confidence levels as in Fig. |l| (referred to the absolute minimum in each panel). 

Concerning the Xp U n minima in Fig. ||, their positions are not particularly interesting for the CI and Ga cases, 
where they are essentially degenerate. More interesting is the case of the SK experiment alone, where the best fit 
(Xpuii = 38.4) is reached at maximal mixing and for Sm 2 — 6.5 x 10~ 10 eV 2 (in agreement with the results in (54|]), very 
close to the QVO coordinates in Table H. Concerning the fit to SNO data only, we find the best fit at Sm 2 — 3.7 x 10~ 5 
eV 2 and tan 2 612 = 0.47, close to the LMA coordinates in Table EI with Xp U n = 25.7. The latter value appears to be 



11 We do not find acceptable solutions in the octant-symmetric vacuum oscillation (VO) regime. For the QVO solution in Table || only 
the highest energy v flux components ($a and •I'hep) have reached the VO regime, while the lowest energy ones are still affected by 
octant-asymmetric quasivacuum effects [ jtq , |49| in the Sun. 

12 Us eful .discus sions of the applications and differences between the goodness-of-fit test and the parameter estimation test can be found 

inggli- r-, 

13 The closest agreement is reached with the global allowed LMA region in Ref. 181. 

14 CHOOZ data are not included in the present analysis, in order to show more clearly the strength of the upper bound on 8m 2 placed by 
solar neutrino data alone. 
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on the lower side of the ±1<t expected range for the x 2 ull m SNO (32 ± a/64)- 15 

Concerning the shapes of the allowed regions in Fig. 0, we note the following facts. None of the experiments 
excludes maximal mixing and Sm 2 — > oo at 99% C.L. This feature is rather well known for Ga and SK data, but 
does not appear in all recent analyses for the CI and SNO data. For instance, the CI contours in Rcf. || appear 
to be more restrictive than ours, which might be due to different estimates of the CI errors. We also note some 
differences between our SNO bounds in Fig. [2] and the SNO official analysis in Fig. 4a of ||: (a) in the QVO region, 
our contours are smooth (as they should); (b) we do not find Q(VO) solutions at maximal mixing for 8m 2 ~ 10 -10 
eV 2 ; (c) our bounds in the LMA region allow maximal mixing and Sm 2 —> oo well within the 99.73% C.L. Concerning 
the points (a) and (b), we think that the differences might depend in J9| on the (numerically delicate) averaging of 
the oscillating terms in the energy or time domain. Concerning the points (b) and (c), some differences might also 
be due to the fact that the analysis in Fig. 4a of || is done without SSM input. Concerning the point (c), we note 
that the exclusion of 8m 2 — ► oo and of tan 2 6*12 = 1 at the 99.73% level by the SNO data alone (as found in Fig. 4a 
of j9|) would be equivalent to the exclusion of the constant P ee — 1/2 case (or of equal 24 and T fluxes, <& e = $^,-7-) 
at the same confidence level. Although the SNO data clearly prefer P ee ~ 1/3 || (see also App. jl|), a 3a rejection of 
P ee ~ 1/2 might be premature. Indeed, from Fig. 3 of |^) it appears that the $ e = Q^.t line touches the 95% error 
ellipse determined by the total SNO rates. Given that the exclusion of relatively large values of 8m 2 and tan 2 9 has 
profound implications in lepton physics (both phenomenologically and theoretically) , we think that the impact of the 
SNO data on such values warrants further investigations. 

From Fig. || it also appears that all experiments largely agree in the LMA region, and that a few QVO "islands" 
below 10 -9 eV 2 also happen to be consistent with all experiments. Moreover, all the experiments appear to be 
generically consistent with some "LOW" or "SMA" regions at least at 99% C.L. However, such regions are somewhat 
different for the different panels in Fig. ^. Concerning the LOW case, there is a reasonable overlap of the CI, Ga, and 
SK bounds at 8m 2 < 10 -7 eV 2 , while SNO prefers 8m 2 > 1CP 7 eV 2 , where the LOW parameters more easily adapt 
to the preferred value P ee ~ 1/3. This tension generates an overall decrease of the LOW likelihood. Concerning the 
SMA case^ the low-tan 2 812 regions separately allowed by CI, Ga, SK, and SNO at the 99.73% C.L. do not overlap. 

Figure |3| allows to appreciate the impact of each experiment in the global fit, by removing one experiment at a time. 
In a sense, Fig. I is the "differ ence" between Fig. [j] and Fig. 0. It can be seen that the removal of either the CI or 
the Ga experiments (upper panels in Fig. ||) weakens the bounds on large mixing, on large values of Sm 2 , and on the 
LOW parameters, but does not alter the situation for vacuum oscillations, which are excluded in both cases. 16 The 
lower left panel in Fig. |^ shows that the removal of the SK experiment would diminish the likelihood of the solutions 
characterized by no or mild energy spectrum distortion (e.g., 8m 2 — > 00 or QVO cases), and, conversely, would make 
some strongly energy-dependent solutions marginally reappear (at small mixing and in the vacuum regime) . Finally, 
the lower right panel in Fig. || shows the pre-SNO situation (but with updated SK, CI, and Ga data), with all the 
well-known multiple solutions to the solar neutrino problem. 

The comparison of Figs. [l]-|3| shows the dramatic impact of SK and SNO in determining the preference for large 
mixing (and especially for the LMA solutions) and the rejection of the SMA solution. However, the CI and Ga data 
still play an important role in determining the shape of the LMA contours, as well as the likelihood of the less favored 
solutions (LOW and QVO), which cannot be rejected on the basis of the present global information. 



C. Separating and grouping pulls 

We discuss the decomposition of Xp U n separate and grouped pulls of observables {x n } and of correlated 
systematics {£ fe } [see Eqs. (|Ls|)~(p2"|)] for the various solutions. 

Table O shows the pulls of each of the 81 observables used in our analysis, corresponding to the four solutions in 
Table 0. 1 ' For the best-fit LMA solution, we also give the range spanned by each pull within the 99.73% C.L. LMA 
region shown in Fig. |l|. It can be seen that, in any case, there are no anomalously large (say, > 3 standard deviation) 
pulls — not even in the SMA solution or at the borders of the LMA region. This fact confirms the results of the 



15 We think that this feature might be partly due to non-optimal binning of the SNO spectrum. Although, at low energy, relatively dense 
binning is required to enhance the effects of the neutral current component, at high energies it is preferable to enlarge the bin width so 
as to match the SNO energy resolution width (analogously to the current SK energy spectrum binning). 

16 In Ref. feci, the Ga impact on the LOW solution has been studied by lowering the total rate from 70.8 ± 4.4 to 66.1 ± 5.3 SNU. The 
authors of ppl] find a corresponding reduction of the Ax^ ow from 6.9 to 3.0, with respect to the LMA minimum. By repeating the 
same exercise, we find a smaller effect: our Ax^qw d ecreases from 10.4 (see Table to 9.0. 

17 The SK and SNO spectrum bins are identified by their energy range and by their nadir interval (day and night bins). See also Apps. O 
and O for notation. 
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TABLE II: Pulls of the obscrvablcs in the various solutions. 
n Observable LMA LMA |min, max | at 3a LOW QVO SMA 



1 


CI rate 










1 — Wa 

— 1.74 


— 2.56, 


— 1 


1 ( 


1 — oc 

— 1.86 


-2.68 


0.14 


2 


Ga rate 










0.39 


— 1.49, 


+ 1 


90 


1.27 





84 


— 0. 19 


3 


G&W - S 










— 1.23 


— 1.23, 


— 1 


22 


— 1.75 


— 1 


25 


— 1.22 


4 


SK 


5 


0,5.5 










— 0.09 


— 0.46, 


+0 


10 


0.14 


— 


01 


0.58 


5 


SK 


5 


5,6.5 




day 


— 0.36 


— 1.06, 


+0 


40 


— 0.06 


— 


52 


0.17 


6 


SK 


5 


5,6.5 




Ml 


— 0.41 


— 0.63, 


— 


29 


— 0.27 


— 


40 


—0.14 


7 


SK 


5 


5,6.5 




M2 


— 1.77 


— 2.00, 


— 1 


64 


— 1.84 


— 1 


72 


— 1.43 


8 


SK 


5 


5,6.5 




M3 


0.03 


— 0.20, 


+0 


26 


0.05 





19 


0.52 


9 


SK 


5 


5,6.5 




M4 


1.11 


+0.93, 


+ 1 


36 


1.46 


1 


30 


1.64 


10 


SK 


5 


5,6.5 




M5 


0.62 


+0.47, 


+0 


86 


0.71 





80 


1.14 


11 


SK 


5 


5,6.5 




core 


— 0.90 


— 1.04, 


— 


67 


— 0.75 


— 


73 


—0.41 


12 


SK 


6 


5,8.0 




day 


1.20 


+0.48, 


+2 


35 


1.15 





56 


0.72 


13 


SK 


6 


5,8.0 




Ml 


1.97 


+ 1.80, 


+2 


04 


1.99 


1 


84 


1.90 


14 


SK 


6 


5,8.0 




M2 


1.31 


+0.82, 


+ 1 


38 


1.08 


1 


29 


1.37 


15 


SK 


6 


5,8.0 




M3 


— 1.29 


— 1.80, 


— 1 


07 


— 1.39 


— 1 


13 


— 1.02 


16 


SK 


6 


5,8.0 




M4 


0.19 


— 0.23, 


+0 


46 


0.62 





41 


0.55 


17 


SK 


6 


5,8.0 




M5 


— 0.81 


— 1.17, 


— 


54 


— 0.78 


— 


60 


— 0.40 


18 


SK 


6 


5,8.0 




core 


— 1.12 


— 1.41, 


— 


87 


—0.99 


— 


94 


—0.73 


19 


SK 


8 


0,9.5 




day 


— 0.40 


— 0.93, 


+0 


56 


— 0.67 


— 


91 


— 1.29 


20 


SK 


8 


0,9.5 




Ml 


0.17 


+0.05, 


+0 


19 


0.13 





09 


— 0.05 


21 


SK 


8 


0,9.5 




M2 


0.24 


— 0.32, 


+0 


42 


0.06 





31 


0.17 


22 


SK 


8 


0,9.5 




M3 


—0.17 


— 0.73, 


+0 


15 


—0.22 





08 


—0.05 


2:5 


SK 


8 


0,9.5 




M4 


1.40 


+0.92, 


+ 1 


76 


1.75 


1 


70 


1.61 


24 


SK 


8 


0,9.5 




M5 


— 0.26 


— 0.65, 


+0 


09 


— 0.20 





02 


0.03 


25 


SK 


8 


0,9.5 




core 


— 0.51 


— 0.83, 


— 


21 


— 0.39 


— 


28 


—0.23 


26 


SK 


9 


5,11.5 




day 


— 0.67 


— 1. 19, 


+0 


19 


— 1.07 


— 


83 


— 1.70 


27 


SK 


9 


5,11.5 




Ml 


— 0.20 


— 0.36, 


— 


13 


— 0.29 


— 


17 


—0.49 


28 


SK 


9 


5,11.5 




M2 


1.01 


+0.45, 


+ 1 


30 


0.88 


1 


22 


0.90 


29 


SK 


9 


5,11.5 




M3 


— 0.55 


— 1.11, 


— 


If) 


— 0.56 


— 


13 


— 0.49 


30 


SK 


9 


5,11.5 




M4 


0.08 


— 0.43, 


+0 


57 


0.40 





55 


0.26 


31 


SK 


9 


5,11.5 




M5 


0. 17 


— 0.21, 


+0 


63 


0.26 





57 


0.45 


32 


SK 


9 


5,11.5 




core 


0.33 


+0.01, 


+0 


69 


0.43 





62 


0.52 


33 


SK 


11.5,13.5 


day 


0.72 


+0.36, 


+ 1 


31 


0.42 





66 


0.11 


34 


SK 


11.5,13.5 


Ml 


0.67 


+0.57, 


+0 


71 


0.61 





71 


0.52 


35 


SK 


11.5,13.5 


M2 


— 2.21 


— 2.64, 


— 1 


98 


— 2.25 


— 2 


00 


— 2.23 


36 


SK 


11.5,13.5 


M3 


— 1.76 


— 2.11, 


— 1 


43 


— 1.72 


— 1 


39 


— 1.65 


37 


SK 


11.5,13. 


5 


M4 


0.36 


+0.02, 


+0 


69 


0.57 





70 


0.55 


38 


SK 


11.5,13.5 


M5 


— 0.80 


— 1.02, 


— 


50 


— 0.72 


— 


59 


— 0.54 


39 


SK 


11.5,13.5 


core 


0.67 


+0.46, 


+0 


92 


0.76 





76 


0.82 


40 


SK 


13.5,16.0 


day 


0.65 


+0.36, 


+ 1 


02 


0.50 


— 


37 


0.48 


41 


SK 


13.5,16.0 


Ml 


0.59 


+0.52, 


+0 


61 


0.56 





25 


0.57 


42 


SK 


13.5,16.0 


M2 


2.03 


+ 1.88, 


+2 


06 


2.03 


1 


79 


2.08 


43 


SK 


13.5,16.0 


M3 


1.51 


+ 1.39, 


+ 1 


57 


1.54 


1 


26 


1.62 


44 


SK 


13.5,16.0 


M4 


1.31 


+ 1. 18, 


+ 1 


39 


1.42 


1 


05 


1.47 


45 


SK 


13.5,16.0 


M5 


— 0.48 


— 0.57, 


— 


41 


— 0.42 


— 


85 


—0.26 


46 


SK 


13.5,16.0 


core 


— 0.22 


— 0.30, 


— 


15 


— 0. 15 


— 


60 


— 0.08 


47 


SK 


16.0,20.0 






0.34 


+0.06, 


+0 


47 


0.31 


— 


19 


0.35 


48 


SNC 


) 


5.0,5 


E 




day 


— 0.99 


— 1.32, 


— 


28 


— 0.49 


— 


66 


—0.08 


49 


SNO 


5.0,5 


5 




night 


0.23 


— 0. 17, 


+ 1 


01 


0.78 





65 


1.19 


50 


SNO 


5.5,6 


0' 


day 


0.94 


+0.63, 


+ 1 


42 


1.23 


1 


06 


1.33 


51 


SNO 


5.5,6 


0' 


night 


— 0.31 


— 0.65, 


+0 


12 


—0.36 


— 


64 


—0.60 


52 


SNO 


6.0,6 


E 




day 


0.09 


— 0.53, 


+0 


64 


— 0.26 


— 


60 


—0.74 


53 


SNO 


6.0,6 


E 




night 


— 0.72 


— 1.55, 


— 


02 


— 1.27 


— 1 


52 


— 1.88 


54 


SNO 


6.5,7 


0' 


day 


— 0.64 


— 1.46, 


+0 


10 


— 1.24 


— 1 


20 


— 1.82 


55 


SNO 


6.5,7 





night 


1.34 


+0.70, 


+ 1 


95 


0.84 


1 


17 


0.40 


56 


SNO 


7.0,7 


5' 


day 


— 0.66 


— 1.24, 


— 


05 


— 1. 14 


— 


57 


— 1.51 


57 


SNO 


7.0,7 


5' 


night 


0.27 


— 0.13, 


+0 


74 


— 0.09 





38 


—0.31 


58 


SNO 


7.5,8 


0' 


day 


0.41 


+0.15, 


+0 


79 


0.16 





37 


0.06 


59 


SNO 


7.5,8 


0' 


night 


0.39 


+0.25, 


+0 


fiN 


0.24 





13 


0.22 


60 


SNO 


8.0,8 


E 




day 


— 0.20 


— 0.26, 


+0 


04 


— 0.27 


— 


65 


— 0.21 


61 


SNO 


8.0,8 


E 




night 


— 0.69 


— 0.73, 


— 


47 


—0.68 


— 1 


20 


— 0.58 


62 


SNO 


8.5,9 


0' 


day 


— 1.31 


— 1.39, 


— 1 


07 


— 1.25 


— 1 


78 


— 1.10 


63 


SNO 


8.5,9 


0' 


night 


0.43 


+0.35, 


+0 


60 


0.51 





27 


0.60 


64 


SNO 


9.0,9 


E 




day 


— 0.72 


— 0.96, 


— 


39 


— 0.51 


— 


20 


— 0.44 


65 


SNO 


9.0,9 


E 




night 


0.54 


+0.07, 


+ 1 


35 


1.03 


1 


09 


1.71 


66 


SNO 


9.5, 10. ( 




day 


— 0.95 


— 1.59, 


+0 


07 


— 0.29 


— 


19 


0.50 


67 


SNO 


9.5,10.( 


1 


night 


— 0.66 


— 1.13, 


— 


07 


— 0.27 


— 


21 


0.17 


68 


SNO 


10.0,10 


5 


day 


0.76 


+0.44, 


+0 


85 


0.80 





75 


0.84 


69 


SNO 


10.0,10 


5 


night 


2.19 


+ 1.70, 


+2 


43 


1.96 


1 


84 


1.75 


70 


SNO 


10.5,11 





day 


0.60 


-0.09, 


+1 


02 


0.20 





18 


-0.14 


71 


SNO 


10.5,11 





night 


-1.76 


-2.54, 


-1 


22 


-2.24 


-1 


93 


-2.60 


72 


SNO 


11.0,11 


5 


day 


0.78 


+0.15, 


+1 


24 


0.41 





99 


0.17 


73 


SNO 


11.0,11 


5 


night 


-0.69 


-1.25, 


-0 


21 


-0.98 


-0 


21 


-1.13 


74 


SNO 


11.5,12 





day 


0.27 


-0.17. 


+0 


69 


0.10 





69 


0.08 


75 


SNO 


11.5,12 





night 


-0.27 


-0.64, 


+0 


10 


-0.33 


-0 


08 


-0.25 


76 


SNO 


12.0,12 


5 


day 


1.28 


+ 1.01. 


+1 


54 


1.31 


1 


17 


1.45 


77 


SNO 


12.0,12 


5 


night 


0.89 


+0.66. 


+1 


14 


1.00 





59 


1.19 


78 


SNO 


12.5,13 





day 


-0.45 


-0.69, 


-0 


13 


-0.26 


-0 


84 


-0.03 


79 


SNO 


12.5,13 





night 


0.64 


+0.44, 


+0 


92 


0.82 





40 


1.02 


80 


SNO 


13.0,13 


5 


day 


1.37 


+ 1.21, 


+1 


60 


1.54 


1 


33 


1.68 


81 


SNO 


13.0,13 


5 


night 


-1.03 


-1.41. 


-0 


57 


-0.64 


-0 


19 


-0.51 
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TABLE III: Separate experimental contributions to the global Xobs ^ or ^ ne var i° us solutions reported in Table jjj, as obtained by grouping 
squared pulls from Table O. 



n 


Experiment 


LMA 


LOW 


QVO 


SMA 


1 


CI (rate) 


3.0 


3.4 


7.2 


0.0 


2-3 


Ga (rate and W - S) 


1.7 


4.7 


2.3 


1.5 


4-48 


SK (44 bins) 


40.4 


42.9 


37.0 


43.1 


49-81 


SNO (34 bins) 


26.2 


28.7 


28.4 


38.5 


1-81 


All (global Xobs) 


71.3 


79.7 


74.9 


83.1 



previous sections, namely, that the allowance or rejection of the various solutions comes from a "collective" effect of 
several experimental observables rather than by a small subset of them. 

Figure |4| provides a graphical version of the first three columns in Table ||, together with an histogram of the pull 
distribution, compared with a univariate gaussian distribution (conventionally normalized to N = 81) to guide the 
eye. It appears that the pulls of the observables have a reasonably symmetrical and gaussian distribution, confirming 
the goodness of the LMA fit at a deeper level than the global \ 2 values. The relatively large pull for the Chlorine 
datum (—1.74) may be regarded as a statistical fluctuation among the others, at the LMA best fit. Notice that the 
CI pull in the SMA solution is instead very small (0.14). Retrospectively, this small pull, together with theoretical 
prejudices against large mixing, appears to be at the origin of a very long detour towards small-mixing oscillations 
in matter. The distributions of pulls of the observables for the QVO and LOW cases (not graphically shown) are 
also reasonably gaussian as for the LMA solution, although with a slightly larger area (given by the corresponding 
Xobs va hies in Table ||). The distribution of pulls for the SMA case (not shown), besides having an even larger area 
(Xobs = 83.1), appears also to be slightly skewed, with an excess of positive pulls. This adds to the statistical problems 
of this solution. 

It is interesting to group the separate experimental contributions to the global Xobs m Table ||, by summing up 



III. show which 



the corresponding squared pulls from Table |I|. The results of this exercise, as reported in Table 
experiment "wins" or "looses" in the various global solutions. In particular, the radiochemical experiments clearly 
win in the SMA and loose in the QVO and LOW solutions, while the fit to the SK spectrum observables appears 
to be rather stable, with only Ax^ bs ~ ±3 variations in the various solutions. The SMA solution tries to make a 
compromise between SK and SNO data, in which SK dominates (having smaller spectral errors), leaving SNO with a 
worse fit (x^bs = 38.5) as compared with the LMA case (Xo bs = 26.2). 

Let us now consider the contributions of correlated systematics shifts to the global fit. Table IV shows the contribu- 
tions of the pulls of systematics in the various solutions. As in Table ^|, for the LMA solution we also show the range 



TABLE IV: Pulls of the systematics in the various solutions. 



k Systematic LMA LMA |min. max| at 3a LOW QVO SMA 



1 


Su 






-0 


05 


-0 


31. 


+0 


42 





19 


0.35 


0.34 


2 


S33 









00 


-0 


13, 


+0 


32 





10 


0.24 


0.18 


3 


S34 









01 


-0 


98, 


+0 


42 


-0 


32 


-0.72 


-0.59 


4 


Si, 14 




-0 


15 


-0 


49. 


-0 


02 


-0 


08 


-0.21 


-0.01 


5 


S17 









38 


-0 


71, 


+ 1 


15 


-0 


35 


-0.47 


-0.83 


6 


Luminosity 







04 


-0 


31, 


+0 


17 


-0 


08 


-0.20 


-0.22 


7 


Z/X 







03 


-0 


89. 


+0 


48 


-0 


35 


-0.68 


-0.60 


8 


Age 









00 


-0 


00, 


+0 


04 


-0 


02 


-0.05 


-0.04 


9 


Opacity 




-0 


05 


-0 


30. 


+0 


52 





22 


0.42 


0.41 


10 


Diffusion 




-0 


02 


-0 


20, 


+0 


43 





18 


0.34 


0.31 


11 


Cbc 






-0 


07 


-0 


22. 


+0 


13 





07 


0.09 


0.16 


12 


iShep 




-0 


03 


-0 


04, 


+0 


00 


-0 


02 


-0.10 


0.02 


13 


8 B v shape 







17 


-0 


71, 


+ 1 


24 


-0 


66 


-0.80 


-1.56 


14 


SK scale 







78 


+0 


51, 


+ 1 


82 





49 


0.49 


-0.31 


15 


SK resolution 





61 


+0 


54, 


+0 


87 





61 


0.06 


0.73 


10 


SK offset 







44 


+ 


33. 


+0 


70 





57 


0.68 


0.34 


17 


SK 


5.0,5.5 




-0 


03 


-0 


18. 


+0 


05 





06 


0.00 


0.27 


18 


SK 


5.5,6.5 




-0 


26 


-0 


61, 


-0 


13 


-0 


10 


-0.28 


0.34 


19 


SK 


6.5,8.0 







54 


+ 


32. 


+0 


67 





70 


0.52 


0.89 


20 


SK 


8.0,9.5 







01 


-0 


00, 


+0 


19 


-0 


08 


-0.03 


-0.42 


21 


SK 


9.5,11.5] 


-0 


14 


-0 


25. 


+0 


21 


-0 


30 


0.14 


-0.76 


22 


SK 


11.5,13.5 


-0 


21 


-0 


31, 


-0 


06 


-0 


29 


-0.10 


-0.45 


23 


SK 


13.5,16.0 





26 


+0 


23. 


+0 


34 





32 


0.11 


0.44 


24 


SK 


16.0,20.0 





01 


+0 


00, 


+0 


02 





02 


-0.01 


0.02 


25 


SNC 


) scale 




-0 


15 


-0 


90. 


+0 


58 


-0 


86 


-1.49 


-1.48 


26 


SNO resolution — 


32 


-0 


41, 


-0 


05 


-0 


16 


-0.52 


0.47 


27 


SNO vertex 





13 


-0 


60, 


+0 


65 


-0 


52 


-0.20 


-1.42 


28 


SNO n capture — 


10 


-0 


40, 


+0 


60 





42 


0.34 


0.94 


29 


SNO n background — 


06 


-0 


27. 


+0 


35 





25 


0.20 


0.55 


30 


SNO LE background -0 


16 


-0 


49, 


+0 


53 





33 


0.28 


0.87 


31 


SNO cross 


section 


04 


-0 


16. 


+0 


21 


-0 


16 


-0.02 


-0.52 
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TABLE V: Separate 


contributions to the global Xs ys 


for the various 


solutions reported in Table h 

L 


. as obtained by 


grouping squared pulls 


from Table |iv[ 










k 


Systematic sources 


LMA 


LOW 


QVO 


SMA 


1-13 


SSM and 8 B u shape 


0.2 


1.0 


2.4 


4.3 


14-24 


SK systematics 


1.7 


1.7 


1.1 


2.9 


25-31 


SNO systematics 


0.2 


1.4 


2.8 


6.6 


1-31 


All (global xt ye ) 


2.1 


4.1 


6.3 


13.8 



spanned by the pulls within the 99.73% C.L. LMA region in Fig. |l|. It appears that such pulls are generally rather 
small, and typically assume minimal values in the LMA solution. In particular, Fig. [5] shows the LMA pull diagram 
for the correlated systematics, none of which presents a significant offset. The smallness of all SSM input offsets (with 
the possible exception of Sn) is particularly impressive. In general, at the LMA best-fit point there is no need to 
stretch any correlated systematics to fit the data. In the LOW and QVO solutions the situation is less ideal (several 
offsets in Table are at the level of ~ 0.5a) but certainly still acceptable. In the SMA case, however, several SNO 
systematic offsets in Table |y| are at the level of ~ la, and the 8 B v shape uncertainty is stretched by ~ 1.5c. Such 
offsets, which act in the direction of reducing the spectral differences between theory and data, produce a significant 
contribution to the SMA fit (xi ys = 13.8). It is quite unlikely that future, possible recalibrations of systematics may 
just happen to cancel out all these offsets, thus giving more chances to the SMA solution. 

It is interesting to group some x 2 contributions of systematics, according to their origin. This exercise is done 
in Table [v|. The SK contribution to \sy S i s rather stable (and small), with a minimum at the QVO solution. This 
situation parallels the SK contribution to Xobs ( see Table |l| and related comments) , and show the pivoting role of 
the SK spectrum in determining the likelihood of all solutions. Conversely, the contributions to Xsy S from SNO and 
from the standard neutrino flux input (SSM and 8 B v shape systematics) increase significantly when passing from 
the LMA to the SMA solution. 



The grand total of the various contributions to the Xp U n fr° m both observables and systematics is shown in Table VI 



This Table shows quantitatively that the LMA solution is in very good agreement with both the experimental data 
and with the SSM. The LOW and QVO solutions provide a slightly less good agreement with the SSM and with 
SK+SNO, and are somehow "borderline" from the point of view of radiochemical experiments. However, they cannot 
be really excluded by any data at present. The SMA solution is instead safely ruled out, mainly as a consequence of 
the SNO data fit (x| NO = 45.1). 





TABLE VI: 


Separate contributions to the 


global Xp ull for the 


various solutions reported in 


Table Q. 




n 


k 


Contributions 


LMA 


LOW 


QVO 


SMA 




1-13 


SSM and 8 B v shape 


0.2 


1.0 


2.4 


4.3 


1 




CI experiment 


3.0 


3.4 


7.2 


0.0 


2-3 




Ga experiments 


1.7 


4.7 


2.3 


1.5 


4-48 


14-24 


SK experiment 


42.1 


44.6 


38.1 


46.0 


49-81 


25-31 


SNO experiment 


26.4 


30.1 


31.2 


45.1 


1-81 


1-31 


All (global x 2 ) 


73.4 


83.8 


81.2 


96.9 



D. Implications for the SSM neutrino fluxes 

In the previous section, we have seen that the LMA solution does not require any significant offset £ fc in the SSM 
input systematics (k = 1, . . . , 12), and that such offsets are relatively small also in the other solutions. Equation [Tt] 
allows to translate such offsets into preferred shift of the neutrino fluxes from their SSM central values. The results 



are shown in Table VII for the various solutions. In the LMA best-fit point, all solar neutrino fluxes are basically 
confirmed: just for illustration, an approximate translation of the LMA flux shifts into variations of the "effective 
solar core temperature" (At c /t c ^ g ||) would formally provide \At c /t c \ < 0.2%. The QVO and SMA solutions 
would formally require —Atc/tc ~ 0.5-1% (a slightly cooler Sun), the LOW case being intermediate between the 
latter and the LMA one. In all solutions, the preferred <&p P values are within a percent from the SSM. The preferred 
(negative) variations of the $b values from its SSM value are instead significant in non-LMA solutions, and appear 
to be consistent with those derived in some $s-free analyses O, O, p3 . This might seem surprising, since we do not 
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TABLE VII: Fractional neutrino flux shifts from the SSM central values (A$i/$j x 100), for the various solutions. 





PP 


pep 


Be 


B 


hep 


N 


O 


LMA 


0.0 


0.0 


0.5 


5.2 


-0.8 


-1.0 


-1.2 


LOW 


0.5 


0.8 


-5.5 


-12.2 


0.0 


-8.3 


-8.3 


QVO 


1.1 


1.6 


-11.5 


-22.2 


-1.8 


-15.1 


-17.0 


SMA 


0.9 


1.3 


-9.9 


-24.2 


+1.8 


-11.4 


-12.8 



exclude any SSM input in our analysis. The reason for such agreement is that the current SNO data implicitly provide 
an experimental determination of $b which is already significantly more accurate than the SSM estimate; therefore, 
there is little difference in making the analysis with or without SSM input for Notice that, in addition, we can 
also quantify the preferred variations of <I>Bo, &n, and $o, which appear to be all negative in non-LMA solutions. 18 
Concerning ^hcp, we do not find any significant preferred variation with respect to Eq. (|l5|). This fact is mainly 
due to the poor sensitivity of the data to this flux (confined to the last few bins in the SK and SNO spectrum) and 
possibly to our careful treatment of SK and SNO spectral uncertainties (see Apps. |C] and [S]). 

Finally, Fig. || shows isolines of the preferred neutrino flux shifts within the LMA solution (superposed at 99% 
C.L.). The variations of $ pp and of 'I'pcp are limited to about a percent, but can be an order of magnitude larger 
for the other neutrino fluxes. In the upper or right part of the LMA region, the current global fits seem to prefer 
variations of the fluxes roughly corresponding to a slightly "cooler" Sun, while in the lower left corner of the LMA 
region the trend is opposite (slightly "hotter" Sun). Should the best-fit LMA value be confirmed by the KamLand 
experiment 13|, the current SSM input would be just "perfect." 



SUMMARY AND CONCLUSIONS 



We have performed a global analysis of solar neutrino oscillations within the 2v active scenario, 19 including all the 
current solar neutrino data and all relevant sources of uncertainties. The statistical analysis has been performed in 
a way which clearly traces the residual contribution of each observable and of each source of correlated systematics 
in the \ 2 fit ("pull" method). It turns out that there is still a multiplicity of acceptable solutions at large mixing 
in the so-called LMA, LOW, and QVO regions (Fig. |l]and Table [|, none of which contradict any data (Figs. || and 
^), as also confirmed by a detailed pull analysis (Tables |l]-|Vl]). In particular, the best-fit LMA solution appears to 
be in very good agreement with all the data (Fig. ||), without requiring any stretching of the correlated systematics 
(Fig. D), contrary to the strongly disfa vored solution at small mixing angle. The striking LMA agreement with all 
standard solar model fluxes (Table VII) is only slightly worsened when moving away from the best fit (Fig. ||). The 
statistical techniques and the treatment of the data underlying these results have been discussed in detail in various 
appendices. 
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APPENDIX A: PROOF OF X 



2 

covar 



— Xpull 



Let us remind that basic ingredients of any \ 2 statistics are: the experimental and theoretical values (i?^ xpt and 
i?^ hcor ) of the N observables to be fitted; the associated uncorrelated errors u„; and the associated set of fully 



18 This information is relevant for prospective studies in the BOREXINO experiment 

19 Three-flavor results will be shown elsewhere. 
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correlated errors c„, due to K independent sources of systematics (k = 1, . . . , K). 

In order to simplify the notation, we normalize both the differences i?° xpt — i?^ hco and the correlated errors to 



the UnS, by defining 



jpexpt Dthe 



(Al) 



and 



Equation (|5|) reads then 



«S = ^ 

Itr. 



iV 

Xcovar = ^ ' A" [finm + X/fc^n^m] ) 
n,m— 1 



(A2) 



(A3) 



while Eq. (m) reads 



Xpuii = mm 



JV 



A" 



*:=i 



where the £fc are gaussian random variables with = and (^) = 1. 



The minimization in Eq. (A4) leads to a set of K linear equations in the unknowns 



(A4) 



K 



N 



h=l 



whose solution is 



K 



N 



£fc = £ Skh £ A " 9 

ft=l n=l 



n i 



where S is the inverse matrix 20 



It turns out that the matrix S is also related to the inversion of the covariance matrix, 

K 

[<5«m + Y^k^nlmT 1 = $nm - £ ^kh q ™ 9 ™ ' 

k,h=l 



(A5) 



(A6) 



(A7) 



(A8) 



Indeed, the product of the above two matrices gives the unit matrix. The above equation reduces the inver sion of the 
N x N covariance matrix in Xcovar to the inversion of the (generally much smaller) K x K matrix in Eq. ([AT]) . 
The last step is to write the Xp U u m t erms of the £ fc 's, 



K 



A p 



a 

k ) 



(A9) 



n=l 



k=l 



and to substitute Eq. (A.6) in Eq. (A9). Expanding the r.h.s. of Eq. (A9), and making use of Eq. (A8), one recovers 
the r.h.s of Eq. (|A3|), namely, 



2 _ 2 
Xcovar — Xpull 



(A10) 



20 Here we deal only with symmetric matrices. Therefore, given a matrix equation such as A = B 1 , we can conventionally write it as 
A nm = [-Bum] 1 without index ambiguity. 
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Finally, we observe that setting 

K 



A„-5>£4, (All) 



•En — 

fe=l 



as in Eq. (g), one gets from Eq. (|A9[) a "diagonal" form for Xp U m 

N K 

xIuk = Ys*1 + YaI > ( A12 ) 

n=l fc=l 

as anticipated in Eq. (|l0|). 

A different proof of the previous relations have been discussed in the context of parton distribution fitting p5| , 
where the pull method is now routinely used |^6|, p7f . 

APPENDIX B: WINTER-SUMMER ASYMMETRY IN GALLEX/GNO 

Earth matter effects can generate an observable winter-summer difference (Rw ~ Rs) m the event rates measured 
in gallium experiments. Such difference can be as large as ~ 6 SNU around the LOW solution p9| |. 
The GALLEX/GNO Collaboration has recently reported the measurement j| (see also [ j55| ) 

R w -R s = -11 ± 9 SNU (GALLEX/GNO) , (Bl) 

whose uncertainty is almost entirely statistical f55f , all systematics being largely cancelled in the difference [Q. In 
our "pull" approach, we simply attach to Rw — Rs an uncorrelated error uws = 9 SNU, with no systematics. 

The definition of "Winter" and "Summer" adopted in (J, [55) is slightly different from the astronomical one, 21 and 
corresponds to 

"Winter" = perihelion (5 jan.) ± 3 months , (B2) 
"Summer" = aphelion (5 july) ± 3 months . (B3) 

Therefore, the solar exposure functions for the above periods, as shown in Fig. fj] in terms of the nadir angle (77), are 
slightly different from the ones given in In particular, the comparison of Fig. (7) in this paper with Fig. 5 in |39[| 
shows the the above definitions lead to a lower (higher) exposure of the innermost trajectories in the mantle during 
winter (summer). The total annual e xpos ure is, of course, unaltered. 



The definitions in Eqs. (B2) and (B3) are particularly useful to smoothly extend the theoretical calculation of 
Rw — Rs from the matter-enhanced regime | p9| down to the vacuum oscillation regime, where seasonal variations 
are instead induced by the eccentricity variations of t he orbital dist ance L (purely geometrical (1/L 2 ) effects being 



factored out). Indeed, with the definitions in Eqs. (|B2j) and (B3), the winter-summer rate difference in vacuum 
happens to coincide with the near- far rate difference (Rn — Rf) previously defined in ]57j . We have then matched 
the results found in |39| and in [ E7J , by considering both matter-induced and eccentricity-induced contributions to 
Rw — Rsi so as to calculate this quantity in the whole oscillation parameter space. The match between the matter and 
vacuum regimes occurs in the quasivacuum range, and is made easier by the lucky circumstance that Earth matter 
effects vanish just when the oscillating terms in the v e survival probability start to be important p8[ . 

The datum in Eq. ( |Bl| ) is compatible with no seasonal asymmetry, adds a slight penalty to the region roughly 
corresponding to the LOW solution (where < Rw — Rs ~ 6 SNU), and modulates the likelihood of the solutions in 
the Q(VO) regime, where both positive and negative values of the asymmetry can occur (—26 < Rw~Rs 1^ +26 SNU), 
the negative ones being slightly favored by the GALLEX/GNO datum. We remark that (positive or negative) seasonal 
effects in gallium experiments, being largely driven by low-energy solar neutrino components, can be compatible with 
the nonobservation of ( 8 B-driven) seasonal effects in SK [[34). 

Finally, we mention that the SAGE experiment has recently reported v event rates in gallium, grouped in 
(bi)monthly intervals ||, which appear to be consistent with no seasonal variations (although no explicit Rw — Rs 
estimate is given in ||). Taking the data in |2| at face value, we argue a slightly positive value for Rw — Rs m SAGE, 



21 Winter and summer were denned in |BS|] as six-months intervals separated by equinoxes and centered, respectively, at the winter solstice 
(21 dec.) and at the summer solstice (21 jun). 
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TABLE VIII: Characteristics of the SK energy spectrum. The first two columns define the bins in terms of observed total electron energy 
E. The third column is the ratio of the SK observed event rates to the SSM jllj predictions, as taken from g, with a slight correction 
to account for our default hep v flux [we take <t>h cp =8.3 rather than 9.3 in units of xlO 3 cm -2 s — 1 , see Eq. (|15[)]. The boron flux is 
<3>5 = 5.05 X 10 6 cm -2 s^ 1 |4l| , as in [||. The fourth column gives the fractional uncorrelated error in each bin, whose only component is 
the statistical error (taken from [Q). The fifth column represents those fractional systematic errors which are uncorrelated among energy 
bins, but correlated in nadir bins (from Q). The sixth, seventh and eighth columns represent our evaluation of the total, boron, and hep 
v event rate in SK (in the absence of oscillation), for 1 kton year (kty) exposure. The ninth, tenth, and eleventh column represent our 
evaluation of the three main systematic effects (given as lcr fractional contributions to R^^) which are fully correlated in energy and in 
nadir, and are generated by uncertainties in the 8 B v energy shape, SK energy scale, and SK resolution width. The last two uncertainties 
also affect the hep contribution to the total rate, as reported in the last two columns (our evaluation, given as la fractional contributions 
to -R^^p)- Finally, we attach an overall systematic offset c off ~ 2.75% ^ to all binned rates R^ co , with full correlation in energy and 
nadir. In the presence of oscillations, the last eight columns (properly separated into nadir bins) are recalculated for each (<5m 2 ,tan 2 812) 
point. 



Bin 


E range 


flm pt /-Rm eor 


u m /R'£ pt 






pthcor 


pthoor 
n m,hep 


100 x c, 




100 x c„ 


/ pthcor 
i,hcp / ^m^ep 


m 


(MeV) 


(B+hop) 


X100 


X100 


(1/kty) 


(1/kty) 


(1/kty) 


8 B shape 


scale 


resol. 


scale 


resol. 


1 


[5.0,5.5] 


0.4672 


8.65 


3.24 


81.795 


81.624 


0.171 


0.52 


-0.07 


-0.19 


-0.38 


-0.09 


2 


[5.5,6.5] 


0.4581 


3.08 


1.43 


140.974 


140.654 


0.320 


0.69 


0.14 


-0.20 


-0.29 


-0.11 


3 


[6.5,8.0] 


0.4730 


1.78 


1.37 


155.433 


155.014 


0.419 


1.05 


0.58 


-0.20 


-0.09 


-0.15 


4 


[8.0,9.5] 


0.4601 


2.02 


1.37 


95.577 


95.239 


0.338 


1.61 


1.32 


-0.09 


0.24 


-0.18 


5 


[9.5, 11.5] 


0.4630 


2.23 


1.37 


60.011 


59.688 


0.323 


2.48 


2.50 


0.00 


0.78 


-0.19 


6 


[11.5, 13.5] 


0.4626 


3.64 


1.37 


18.197 


18.004 


0.193 


3.87 


4.50 


1.59 


1.67 


-0.07 


7 


[13.5, 16.0] 


0.5683 


6.88 


1.37 


3.768 


3.660 


0.108 


5.82 


7.47 


4.51 


3.14 


0.42 


8 


[16.5,20.0] 


0.5637 


26.3 


1.37 


0.245 


0.212 


0.033 


8.09 


12.32 


11.26 


6.09 


2.41 



a The total event rate /£ thco is here normalized to the efficiency-corrected, no-oscillation value of 556 events/kty, as graphically derived from [ 
with a correction for the updated SSM $B,hcp values. This specific value for eo is unimportant in practice, as far as the ratios in the 
column are used in the analysis. 
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with a total uncertainty comparable to that of GALLEX/GNO. The combination of a slightly positive (SAGE) and 
a slightly negative (GALLEX/GNO) winter-summer difference would then provide a central value closer to zero for 
R\v — Rs, very consistent with the LMA solution J39|. An official evaluation of Rw — Rs from SAGE [to be combined 
with the one in Eq. (Bl) from GALLEX/GNO] appears thus a desirable input for future analyses. 



APPENDIX C: THE SK ENERGY-NADIR DIFFERENTIAL SPECTRUM 



The SK energy-nadir spectrum of electrons induced by neutrino elastic scattering is fundamental to constrain the 
solar neutrino parameter space. Therefore, we think it useful to describe in some detail our improved error estimates 
and theoretical calculations. 

Concerning the "error format" , the pull method used in this paper does not leave any freedom in the % 2 treatment 
of the spectrum, which is uniquely defined by providing, for each bin rate R n , the uncorrelated error component u n , 
and the fully correlated error components due to independent fc-th sources of systematic errors. 22 The task is thus 
reduced to a careful evaluation of such components. To reach this goal, we combine information from SK j^] and from 
our own evaluation of systematics. 

Table VIII shows the main characteristics of the SK binned energy spectrum. The energy bins 2-7 are further 
divided into seven nadir angle bins |(| , reported in Table IX . In each energy-nadir bin, the statistical error represents 
the only SK uncorrelated error component (u n ). All other error sources in SK are correlated among bins, at least 
in nadir. Indeed, apart from the obvious SSM systematics, the SK bin rates are further affected by 8 + 3 + 1 = 12 

correlated systematics, as discussed in the following. 

The first eight systematic errors, listed in the fifth column of Table VIII, are uncorrelated among energy bins. In 



each energy bin, however, they are fully correlated in nadir [^) . Such errors represent those (nadir-independent) data 
reduction uncertainties which are specific of each energy bin, independently of the others. 

Three further systematics, fully correlated among all energy-nadir bins, are induced by 8 B v spectrum shape 



22 Therefore, we abandon our previous \ 2 approach in terms of separated SK total rate and spectral shape information |58 
correct, cannot be exactly cast in a "pull" form. 



which, although 



16 



TABLE IX: Fractional contributions to the total event rate (R/Rtot, our evaluation), in the case of no oscillation, from each of the seven 
SK nadir intervals, identified as "day" bin, "mantle" bins (M1-M5), and "core" bin. The observed SK event rates in each energy-nadir 
bin, together with their statistical errors, are reported in 

Nadir bin cos rj range i?/i?,t Q t 

Day [-1,0] 0.5000 

Ml [0,0.16] 0.0685 

M2 [0.16,0.33] 0.0777 

M3 [0.33,0.50] 0.0984 

M4 [0.50,0.67] 0.1025 

M5 [0.67,0.84] 0.0839 

Core [0.84, 1] 0.0690 

Day+Night [-1-1] 1.0000 



uncertainties jl6| and by the SK energy scale and resolution uncertainties ||. We note two facts. Firstly, these 
sources of systematics act differently on the B and hep components in SK: for instance, the hep component is obviously 
unaffected by the 8 B shape uncertainty. 23 Secondly, the relative error signs (for each source of systematics) are relevant: 
for instance, an increase in the energy resolution width flattens the — steeply falling — SK energy spectrum, the high- 
energy (low-energy) part being thus enhanced (suppressed), as evident from the eleventh and thirteenth column in 
Table |VIDj . Since the information about the relative error signs and about the separate B and hep components is not 
given in |6j], we perform our own evaluation as follows. 

The 8 B spectrum shape error is evaluated by attaching to the default neutrino spectrum [[l7j the ±la shape 
perturbations evaluated in fl46|| , and calculating the corresponding fractional variations in the absolute rates. 24 The 
energy scale uncertainty is evaluated by shifting the centroid of the (gaussian) energy resolution function by 0.64% 
H , namely, by taking T' — > T"(l + 0.0064), where T" is the true electron kinetic energy (the quantity calibrated 
in SK). The resolution uncertainty is evaluated by perturbing the energy resolution width o~t (see Eq. E3) by 2.5% 
||, <jt or(l + 0.025). In all cases, the same uncertainties are calculated for opposite perturbations, and small 
error asymmetries are averaged. The results are given in the last five columns of Table VIII, where the contributions 
from 4>s and $hop are separated. As far as the B flux component is concerned, there is reasonable agreement in size 
with the corresponding SK error evaluation in except for the 8 B v spectrum shape. We are unable to explain 
such difference. Our error assignment is completed by an overall SK systematic offset (2.75%, symmetrized from ||), 
which is attached to all energy-nadir theoretical rates with full correlation. This error mainly represents the overall 
uncertainty of the data reduction efficiency ]54| , affecting the whole spectrum. 

Notice that, in Table VIII , the last eight columns refer to the no oscillation case. In the presence of oscillations, one 
usually updates only the theoretical rates R^ 00 , and assumes that the fractional errors in the last five columns are 
approximately unchanged. The latter assumption is a very good approximation for spectra with no or mild distortions, 
but is not strictly applicable in the whole oscillation parameter space. Therefore, for t he sake of accuracy, in our 
oscillation analysis we recalculate the fractional errors in the last five columns of Table VIU , for any (Sin 2 , tan 2 #12) grid 
point (and for any nadir bin). This improvement leads to minimal differences in the local Xp U n minima (corresponding 
to almost undistorted spectra), but is not totally negligible in deriving the borders of, say, the 99.73% C.L. allowed 
regions, where it can lead to |A% 2 | ~ 1 differences. 

Concerning the calculations of the v e survival probability, we compute Earth matter effects | f30| t hrough eight 
relevant shells of the Preliminary Reference Earth Model (PREM) plj, namely (in the language of ]61|): (1) Ocean; 
(2,3) Crust layers; (4) LID + low velocity zone; (5) Transition zone; (6) Lower Mantle + D-zone; (7) Outer core; and 
(8) Inner core. In each of them, we approximate the radial density profile through a biquadratic parameterization, 
which allows a fast and accurate analytical calculation of the relevant transition probabilities in the Earth [p0| . 

In the (quasi)vacuum oscillation regime, the oscillating term in the v e survival probability P ee depends implicitly 
upon the daily time through the orbital distance L(Td), where = 2-7T • day/365, with = at winter solstice. 25 As 
noticed in |62f , the different exposure of the "day" and "night" bins in terms of induces a slight day-night vacuum 



Analyses where the &hep/^B ratio is taken as a free parameter should separately rescale the different <E>b and "I>hep error components 
in the evaluation of the total error. In fact, one cannot calculate first the total (<E>s + $hcp) systematic error assuming the SSM ratio for 
^hep/^Si and then use the same error when such ratio is significantly varied (say, by factors ~ 10). We think that this remark should 
be taken into account, when placing upper bounds on "I'hcp from the analysis of the high-energy tail of the SK spectrum. 
Conventionally, we denote as +fcr perturbation of the 8 B v shape the one which moves the energy spectrum to higher energies. 
The purely geometrical 1/L 2 (r,j) variation is assumed to be already corrected for in the SK data. 
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difference in the time-averaged P ee . Concerning the SK energy-nadir spectrum, we take into account the different 
exposures E(r)x, 772) for each nadir bin range [rji , 772] in Table [IX| as follows. In a given day of the year (r^), the 77 range 
spanned by a detector at latitude A (equal to 36.48° for SK) is [77mm (jd), 77 max (Td)], where r\xDxrSjd) = A + <5s(r<j), and 
T/max^d) = 7T — A+(5s(t c ;), having defined the solar declination ftsijd) through sin (5s = — sini cost^, with sini = 0.3978 
(inclination of Earth axis). The daily exposure function E for a generic [771,772] nadir bin is then 

£■(^,771,772) = [T h (T d , min{r72, 77 max }) - T h (r d , max{?7i, ?7min})] /V , (CI) 

where (7^,77) is the hourly time (normalized to 2tt and starting at midnight) corresponding to the angle 77 during 
the day T d , 26 

(COS 77 \ 
; — h tan Atan<5s . (C2) 
cos A cos d s J 



Figure g shows the exposure functions E(r d ) for the seven SK nadir bins listed in Table IX, as well as for the SNO day 
and night bins (the functions sum up to unity in both cases). It can be seen, as intuitively expected, that the SK core 
bin is sensitive only to extreme values of L{r d ) (close to the orbital perihelion and aphelion), while the two outermost 
mantle bins (Ml and M2) are almost equally sensitive to all values of L(T d ), being crossed by solar neutrinos during 
the whole year. We take into account the different bin exposures in Fig. || when time-averaging the oscillating terms 
in the (quasi)vacuum regime, for both SK and SNO. 

The results of our oscillation analysis of the SK spectrum are given in the SK panel of the (previously discussed) 
Fig. ||. The C.L. contours compare well with the official SK ones, as reported in Fig. 2 of M (we also obtain very 
good agreement for the normalization- free case, not shown). 

A final remark is in order. All the SK systematic error sources discussed in || (and adopted in this work) are fully 
correlated in nadir. This means that no allowance is given in [pi for uncorrelated nadir shape variations, apart from 



the obvious statistical fluctuations. However, it is known from previous SK publications [22, |63j and presentations 
p9[ that some systematics do not cancel in day-night differences, and may thus, in general, affect the nadir bins in 
uncorrelated ways. In addition, Earth matter density and chemical composition uncertainties may also contribute to 
small (and different) errors [ )S4"| , |65| in the various nadir bins p(| via v e regeneration effects. In particular, one should 
not forget that the widely used PREM model is spherically symmetric by construction, and that local departures from 



the world-averaged density are to be expected. For instance, the first half of the SK mantle bin Ml (see Table IX) 
is sensitive to the PREM "Ocean+Crust" density, whose local characteristics at the SK site may well be different 
from the world average. The inclusion of such additional uncertainties might give a little more freedom to fit possible 
distortions in the nadir distributions. Therefore, if the (currently weak) hints for an excess of night events in SK [5^| 
and in SNO will be corroborated by future data, an improved discussion of the nadir spectrum uncertainties will 
be useful to precisely assess their statistical significance. 

APPENDIX D: TREATMENT OF THE SNO DAY-NIGHT ENERGY SPECTRUM 

The solar neutrino events observed in the SNO spectrum Q cannot be currently identified as being of ES, CC, 
or NC type on an event-by-event basis. A statistical separation is possible, however, by exploiting their different 
distribution in terms of suitable variables, the most important being the observed (effective) electron kinetic energy 
T. Since the ES and CC distribution shapes in T depend upon the oscillation parameters, also the inferred CC and 
NC rates depend on such parameters, as stressed in [io|. 

In the oscillation analysis, however, it is not necessary to perform a separate fit to the ES, CC, and NC components 
as done in [Q. Given the oscillated predictions for these components, one can simply add them up (together with 
the known background rates) in each energy bin, calculate the corresponding systematics, and fit the observed SNO 
day-night energy spectrum. This method, which has been dubbed "forward fitting" by the SNO Collaboration [f!o|| , 
allows to take into account the full spectral information (central values and errors of each bin). In the following, we 
describe our implementation of such method in the Xp U n evaluation. 

Table [X] reports some relevant characteristics of the SNO spectrum, including our evaluation of the SNO neutrino 
signal components and their errors (small error asymmetries being averaged out). The effects of the 8 B v spectrum 
shape error (not included in the official SNO analysis ^ ||) are estimated as in SK (see App. |^). Notice that a +1<t 



26 Any ang le r) € [r7 m i n , t| maI ] is spanned twice each day. This explains the appearance of n (rather than of 2n) in the denominator of 
Eq. Ef). 
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TABLE X: Characteristics of the SNO spectrum, divided in 17 day and 17 night bins (first two columns). The experimental rates |g are 
given in the third and fourth columns, and include signal (CC+ES+NC) and background events. Our evaluation for the corresponding 
CC, ES, and NC components and of their 8 B v shape, energy scale and resolution errors (in the absence of oscillation, day+night) are 
given in the remaining 12 columns. The NC signal is negligible in the last 8 bins. All the quoted numbers refer to the total + <t>hcp 
contributions. In the analysis, however, the and 'I'hep contributions and their errors are separately evaluated, analogously to the SK 
spectrum. In the presence of oscillations, the last twelve columns are split into day and night contributions and updated in each point of 
the parameter space. Backgrounds are oscillation-independent and are treated according to pcj. 



Bin 

m 


T ranj^c 
(MeV) 


^cxpt 
day 

(l/kty) 


"night 

(l/kty) 


nthcor 

"cc 
(l/kty) 


fractional errors X 100 
shape scale rcsol. 


„thoor 

(l/kty) 


fractional errors X 100 
shape scale rcsol. 


otheor 

(i/kty) 


fractional errors X 100 
shape scale rcsol. 


1 


[5, 5.5] 


351.1 


399.7 


822.2 


0.42 


-2.30 


0.10 


135.8 


0.64 


0.26 


-0.26 


293.8 


1.55 


0.88 


-4.33 


2 


[5.5, 6] 


330.9 


313.4 


875.6 


0.63 


-1.89 


-0.21 


121.7 


0.77 


0.57 


-0.35 


248.5 


1.55 


3.73 


-2.84 


3 


[6, 6.5] 


299.6 


272.2 


903.0 


0.86 


-1.40 


-0.53 


107.8 


0.91 


0.93 


-0.46 


172. 1 


1.55 


7.05 


0.43 


_j_ 


\P. ^ 71 
[0.0, l\ 


29 9 _J 


OA Q fi 


on9 n 


1.12 


-0.85 


-0.83 


94 1 


1.07 


1.32 


-0.54 


Q7 7 


1.55 


10.85 


5.43 


5 


[7,7.5] 


191.2 


235.0 


872.8 


1.41 


-0.22 


-1.09 


81.0 


1.24 


1.77 


-0.58 


45.4 


1.55 


15.16 


12.18 


6 


[7.5,8] 


148.9 


176.6 


817.7 


1.73 


0.51 


-1.30 


68.7 


1.44 


2.28 


-0.58 


17.3 


1.55 


19.96 


20.68 


7 


[8,8.5] 


128.7 


122.1 


740.9 


2.08 


1.33 


-1.42 


57.2 


1.66 


2.85 


-0.54 


5.4 


1.55 


25.26 


31.02 


8 


[8.5,9] 


139.7 


134.1 


648.5 


2.47 


2.25 


-1.41 


46.9 


1.90 


3.48 


-0.43 


1.4 


1.55 


31.19 


43.38 


9 


[9,9.5] 


90.1 


95.6 


547.4 


2.89 


3.29 


-1.24 


37.7 


2.16 


4.19 


-0.23 


0.3 


1.55 


37.69 


58.11 


10 


[9.5, 10] 


82.7 


86.3 


444.9 


3.36 


4.45 


-0.85 


29.6 


2.45 


4.99 


0.07 











11 


[10, 10.5] 


66.2 


62.4 


347.5 


3.86 


5.74 


-0.18 


22.8 


2.77 


5.87 


0.51 











12 


[10.5, 11] 


49.6 


59.8 


260.4 


4.40 


7.16 


0.82 


17.1 


3.12 


6.86 


1.13 











13 


[11, 11.5] 


31.3 


41.2 


186.8 


4.97 


8.70 


2.21 


12.5 


3.51 


7.94 


1.97 











14 


[11.5, 12] 


18.4 


21.2 


128.2 


5.56 


10.37 


4.04 


8.9 


3.92 


9.14 


3.07 











15 


[12, 12.5] 


9.2 


18.6 


84.0 


6.17 


12.15 


6.35 


6.1 


4.36 


10.45 


4.47 











16 


[12.5, 13] 


11.0 


15.9 


52.6 


6.79 


14.01 


9.13 


4.1 


4.83 


11.87 


6.23 











17 


[13,20] 


9.2 


9.3 


71.1 


7.68 


17.83 


16.34 


6.5 


5.97 


15.71 


12.48 












shift of this systematic uncertainty increases the number of signal events, especially at high energies. The fractional 
increase is obviously constant for the NC events, which have no memory of the original v energy. From Table [x| it 
appears that, in the bulk of the SNO spectrum, the 8 B v shape uncertainties are not negligible as compared with 
the (purely instrumental) SNO energy calibration and resolution uncertainties. The latter two errors are evaluated 
by shifting the centroid of the energy resolution function and varying its width in the way described in Ref. J40[ , 
and calculating the corresponding variations in each bin of the CC, NC, and ES spectrum. Actually we separate the 
<&b and $hcp components (not shown in Table |x|), as for the SK spectrum. In the presence of oscillations, we not 
only update the day and night rates in each bin, but also their fractional systematic errors, for each point in the 
mass-mixing plane. As for SK, the effect of this improved estimate of fractional systematic errors is rather small in 
the local % 2 minima, but increases towards the borders of the 99.73% C.L. regions, where the spectral distortions can 
be more sizeable. In the fit to SNO data only, such improvement can lead to Xp un variations as large as \A\ 2 \ ~ 4 
at the 3cr borders, and is thus not totally negligible in deriving precise C.L. contours for the less favored solutions. 
Finally, notice that the three systematic error sources in Table |X| can induce, in the various bins, event rate variations 
with different relative signs, which are taken into account in the analysis. 

Besides the previous three systematic error sources, we include a vertex reconstruction uncertainty, whose +1<7 
effect is to increase the CC and ES rates by +3% and the NC rates by +1.45% Q. The CC and NC binned rates are 
also affected by a systematic cross section uncertainty, whose +la effect is to increase them by +1.8% and +1.3%, 
respectively [||. 27 Finally, the various SNO background components are affected by the so-called neutron capture, 
neutron background, and low-energy (LE) background correlated systematics, whose oscillation-independent effects 
(included in our analysis) are reported in p0| and not repeated here. 

When all such inputs are included, our fit to SNO data only provides the allowed regions given in the lower right 
panel of Fig. |[ Such regions appear to be somewhat different from those found in Ref. H, as discussed in Sec. IV B. 



The correlation between the NC and CC cross section errors can be safely taken ~ 1 in the SNO analysis. 
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APPENDIX E: QUASI-MODEL-INDEPENDENT ANALYSIS OF SK AND SNO 



In this last appendix we present a new version of the model-independent comparison of SK ES and SNO (CC and 
NC) total rates, first proposed in |67j and then applied in |38| to earlier SNO CC data. 

For the current SNO threshold T SNO > 5 MeV [§, the SK threshold equalizing the SNO CC and SK ES response 
functions 



turns out to be Tsk > 6.74 MeV, corresponding to a total electron energy E$k > 7.25 MeV, 
SNO CC response (T SNO > 5 MeV) SK ES response (E SK > 7.25 MeV) 



(El) 



In calculating the response functions, we have taken into account the latest evaluation of the energy resolution width 
a T in SNO M, 



cr T (SNO) = -0.0684 + 0.331^ + 0.0425 T' , 



(E2) 



where T" is the true electron kinetic energy, and both T" and <jt are given in MeV. We have cast, in the same form, 
the latest evaluation of cr T for SK (graphically reported in ) , 



ctt(SK) = 0.25 + 0.20%/T 7 + 0.06 T 1 . 
From the spectral information discussed in App. we estimate (derivation omitted) 

$lf ~ 2.35 ± 0.029 ± 0.080 ± 0.045 (E > 7.25 MeV) , 



(E3) 



(E4) 



in units of 10 6 cm -2 s _1 . In the above equation, the first error is statistical, while the second represents the (properly 
propagated) sum of systematic errors, with the exception of the 8 B shape uncertainty, given separately as a third 
error component. 

While the above SK ES flux estimate is safe for model-independent analyses, the current SNO CC flux || is not, 
since its value does depend on model-dependent assumptions, governing the CC spectrum distortions and thus the 
CC/NC separation Efl ]. However, we can resort to a quasi-model-independent analysis, by assuming that the only 
observable effect of new neutrino physics can be embedded in a shift of the first spectral moment p9| , namely, in a 
linear distortion (tilt) of the CC energy spectrum. 28 This assumption is reasonably general, since the current SNO 
statistics is not high enough to really constrain higher moments, and since we know from SK that only the scenarios 
with mild distortions can survive. In this case, the normalization-preserving form of a generic linear distortion for the 
CC differential energy spectrum reads 



dN CC {T) dNcc (T) 



dT 



dT 



1 



T — (T) 

(t 2 ) - (ry 



A(T) 



(E5) 



where (T) and (T 2 ) are the first and second moment of the undistorted spectrum, and A(T) is the shift in the first 
moment. 

Assuming A(T) free, our fit to the SNO energy spectrum gives the following CC and NC fluxes, 



$£ N C U = 1.78±0.10(stat.) , 
= 4.90±0.66(stat.) , 



<T>SNO 



(E6) 
(E7) 



with correlation p = -0.8. By setting A(T) = 0, we would instead obtain = 1.73 ± 0.06(stat.) and = 

5.29 ± 0.43(stat.) with correlation p = —0.62, in good agreement with the results of || ^| for undistorted spectrum. 
The larger statistical errors in Eqs. (E6) and (E7) are the price to pay to allow possible linear distortions in the fit. 

Systematic errors are attached as follows. From the (positively and negatively correlated) SNO error components 
reported in Table II of we estimate the (CC,NC) experimental systematics as (5.2%, 8.8%), with correlation 
p = —0.5. The corresponding theoretical cross section uncertainties [(1.8%, 1.3%) from Table II of ||] are assumed to 
have p ~ 1. Finally, if we repeat to the fit leading to Eqs. (E6) and (E7) with a +la perturbed 8 B neutrino energy 
spectrum, we obtain the variations (—1.1%, +3.3%), which are fully correlated among themselves and with the third 
SK error component in Eq. (E4). 



The numerical results for a linear distortion of the CC component of the SNO spectrum are very similar to those obtained for a linear 
distortion of the CC+ES components (not shown), since the ES contribution is relatively small in SNO. 
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In conclusion, we get (in units of 10 6 cm 2 s 1 ), 



(T> SK 
^ES 

$SNO 



-cc 

SNO 
NC 



2.35 ±0.10 
1.78 ±0.14 
4.90 ±0.80 



with correlation matrix 




(E8) 
(E9) 
(E10) 



(Ell) 



The above SK and SNO equalized fluxes can be useful to constrain generic models of new physics (alternative to — 
or coexisting with — usual mass-mixing oscillations), as far as their main distortion effect on the SNO CC spectrum 
is approximately linear in T. Within this quasi-model-independent assumption, and in the hypothesis of active 
oscillations, the fluxes in Eqs. (E8)-(E10) are still linked by the exact relations pTJ 



<I> SK - 



"ES 
SNO 
CC 
SNO 
NC 



$B[(Pee)+r(l-(P ee ))} , 
$ B (Pee) , 



(E12) 
(E13) 
(E14) 



0.154 is the ratio of v^ T and v e CC cross sections 



where P ee (E v ) is not necessarily constant in energy, and r 

(averaged over the current SNO equalized response function). In the above equations, $b is the true S B flux from the 
Sun (generally different from the SSM value), and (P ee ) is the energy average of P ee (E v ) over the response function 
(equal in SK and S NO) 

Equations (Eg) -( Ell ) and (E12)-(E14) overconstrain the two parameters $s and (P ee )- Figure || shows, in 
(<&b, (Pee)) coordinates, both separate and combined bounds at the 2a level for e ach datum (A% 2 = 4). There 



is very good agreement between any two out of the three data in Eqs. (|E8|) -( E10 ). Their combination strength- 
ens previous model-independent indications |38| ^8) for a consistency of $b with the SSM prediction and for 
(Pee) ~ 1/3 in the SK-SNO energy range. 
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FIG. 1: Global results of the solar neutrino data analysis, including 81 observables and 31 sources of correlated systematics. 
The parameter space (5m 2 , tan 2 812) refers to the scenario of 2u oscillations among active states. The relevant \ 2 minima in 
the LMA, LOW, and QVO regions are given in Table A 




FIG. 2: Results of the solar neutrino data analysis, as obtained by separating four classes of observables: (i) the chlorine 
rate; (ii) the average SAGE+GALLEX/GNO gallium rate plus the GALLEX/GNO winter-summer difference; (iii) the SK 
energy-nadir spectrum; and (iv) the SNO day-night spectrum. 




FIG. 3: Results of the solar neutrino data analysis, as obtained by excluding each of the four data sets in Fig. 
global set used in Fig. Ill 
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FIG. 4: Diagram of pulls {x„} n =i,...,si for observables at the LMA best-fit point. See the text for details. 
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FIG. 5: Diagram of pulls {£fe}fc=i ... 31 for correlated systematics at the LMA best-fit point. See the text for details. 
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FIG. 6: Preferred deviations of the solar v fluxes from their SSM central values [|4l|], in the region of the LMA solution 
(superposed at 99% C.L.). Dotted (solid) isolines roughly correspond to a slightly "cooler" (hotter) Sun. See the text for 
details. 
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FIG. 7: Solar exposure functions (at the Gran Sasso latitude) for the "winter" and "summer" periods defined in Eqs. (BS) 
and (Si). 
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FIG. 8: Daytime solar exposures of each SK nadir bin and of the SNO day and night bins. These exposure functions are used 
for accurate time averages of the oscillating terms in the (quasi)vacuum regime. 
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FIG. 9: Quasi-model-independent comparison of SK (ES) and SNO (CC and NC) total event rates in the plane charted by the 
boron flux <E>b and by the average Ve survival probability {Pee), for equalized SK and SNO response functions. The evaluation 
of the SNO rates includes possible linear distortions in the energy spectrum. The 2a bands for each datum (A^ 2 = 4) appear 
to be in very good agreement with each other for P ee ~ 1/3 and for <J>b close to its SSM prediction |H|| . The combination of 
the 2a bands is also shown (slanted elliptical region). 



